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It is shown that nearest-neighbor antiferromagnetic interactions of identical Ising spins on imbalanced 
bipartite lattice and imbalanced bipartite hierarchical fractal result in ferrimagnetic transition instead of anti- 
ferromagnetic one. As evidenced by the existing low-field experiments, dilute Ising antiferromagnets may also 
experience such ferrimagnetic transitions due to the imbalanced nature of the magnetic percolation cluster 
when it coexists with the percolation cluster of vacancies. 
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The system of the identical Ising spins on the sites 
of some crystalline lattices with the nearest-neighbor 
antiferromagnetic (AF) exchange may have magnetized 
ground states. In such states there would be antipar- 
allel neighboring spins, as interaction dictates, but the 
whole numbers of up-spins and down-spins would dif- 
fer. One such 2d lattice is shown in Fig. 1(a). Here 
two sublattices with parallel up and down spins in the 
ground state are shown by filled and empty circles cor- 
respondingly. We see that in the unit cell there are 
one filled circle and two empty ones so we get ±1/3 
magnetizations in two globally-reversed ground states 
for nearest-neighbor AF on this lattice. Thus this AF 
model has a couple of ferrimagnetic ground states with 
Iboth staggered L = ({Sa) — (Sb)) /2 and homogeneous 
M = (2 (Sa) + (Sb)) /3 magnetizations. 

One can easily show that this ordering persists up 
to finite T c . Summing the Gibbs function over spins on 
sublattice A (empty circles) we get the Gibbs distribu- 
tion for the spins on the sublattice B having effective 
ferromagnetic Hamiltonian. Indeed, for each link with 
Sa spin we have (J is AF exchange) 

cxp [S A (Sb + S' B J/T)] = 2 cosh [(S B + S' B J/T)] 

Sa=±1 

= 2exp[K B (S b S'b + 1)] , 2Kb = In cosh (2 J/T) 

Hence, the ordering of Sb spins is described by the ferro- 
magnetic Ising model on the square lattice, so (Sb) i= 
for K B > § In (y/2 + l) ffl or 



T <T C = 2 J I In 



As (Sb) is a linear combination of L and M, the ordered 
phase is a ferrimagnetic one. This implies that homoge- 
neous magnetic field H has a part conjugated with the 
order parameter so magnetic susceptibility diverges at 

T L and the transition b ecomes smeared at finite H. 
y e-mail: pntim@live.ru 
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Fig.l. Examples of imbalanced bipartite graphs with 
different numbers of sites in sublattices A (open cir- 
cles) and B (filled circles), Na > Nb- (a) - fragment 
of regular 2d lattice, dotted line shows the unit cell; (b, 
c, d) - clusters of dilute square lattice. In the ground 
state short-range Ising AF on them would have paral- 
lel spins on A and B sublattices and, hence, a nonzero 
magnetization. 

Thus we have a simple example showing that 
nearest-neighbor AF interaction of identical Ising spins 
may result in the macroscopic ferrimagnetic order. This 
is in apparent distinction with conventional ferrimag- 
nets having several different magnetic moments in a cell. 
It may look rather exotic in the realm of real crystals yet 
such situation can be frequent in disordered Ising AF, 
first of all, in dilute Ising AF on bipartite lattices. These 
lattices can be divided in two subsets of sites, A and B, 
such that all bonds are of the A-B type, i. e. there 
are no bonds inside A and B subsets [2]. Apparently, 
the Ising AF on such lattice is non-frustrated having all 
spins up on sublattice A and down on sublattice B or 
vice versa in its two degenerate ground states. Their 
magnetizations are 
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Here Na and Nb are the numbers of sites in A and B 
sublattices and we choose rj < 1 for definiteness. 

Seemingly, all known non-frustrated AF crystals 
with just one sort of magnetic ions have bipartite lat- 
tices that are the balanced ones, that is with r\ = 1 and 
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purely AF ground states, while Fig. 1(a) shows the im- 
balanced bipartite lattice with 77 < 1 (77 = 0.5). Yet 
the dilution of balanced bipartite lattices results in ap- 
pearance of a number of isolated clusters, mostly with 
77 < 1, those with 77 = 1 being the rare exceptions. Figs. 
l(b, c, d) show the imbalanced clusters on the square 
lattice. So at T = and arbitrarily small magnetic field 
dilute bipartite AF must show nonzero magnetization 
due to the presence of such imbalanced finite clusters. 
This circumstance was first noticed by Neel [3J. Still 
it stays unnoticed that for some concentrations of mag- 
netic ions p the giant percolation cluster may also have 
the average imbalance ratio rj p < 1. 

Indeed, in finite sample the role of percolation clus- 
ter belong to that with the largest number of sites and 
very probably it is imbalanced, as most of them. How- 
ever, in the thermodynamic limit (N — > 00) r\ p will tend 
to unity if there are only finite clusters of vacancies. Ap- 
parently, such finite clusters cannot make infinite lattice 
imbalanced as for every cluster deleting unequal number 
of sites from A and B sublattices there exists (with the 
same probability) the shifted cluster of the same form 
which restores the balance. Thus at 1 — p c < p < 1 
the imbalanced percolation cluster can only exist as a 
finite-size effect. Meanwhile, at p c < p < 1 — p c there 
is infinite percolation cluster of vacancies to which this 
argument does not apply. Hence, here rj p < 1 may also 
hold in the N — > 00 limit in some crystal lattices. Then 
the ground state magnetization of dilute AF in this in- 
terval will be 

1 

m p (H = +0) = fi-Hw p (r,)d V + \-^ 

J 1 + V 1 + Vp 



Here W p (77) is the imbalance distribution function of fi- 
nite clusters. Now it seems that neither W p (77) nor ?7 P 
are known for the crystal lattices. So to find them is 
quite relevant task for the physics of dilute short-range 
AF. 

The magnetization of finite clusters vanishes at fi- 
nite temperatures, but that of the percolation cluster 
could persist up to a finite T c making the transition a 
ferrimagnetic one as in the regular imbalanced lattice 
described above. Then macroscopic features of dilute 
AF would be those of ordinary ferromagnet in spite of 
the presence of antiparallel neighboring spins in the or- 
dered phase. In such a case, quite probable for dilute 
AF in a field (DAFF) at p c < p < 1 — p c , the mapping of 
this model onto random-field Ising magnet (RFIM) [?] 
would be no longer valid. Indeed, the derivation of this 
correspondence in Ref. [4] is made under the assump- 
tion that DAFF has a transition into AF phase which 



does not hold if, instead, there is a ferrimagnetic tran- 
sition (smeared at finite H). The evidences in favor of 
DAFF ferrimagnetism can be found in experiments on 
several dilute Ising AF with p < 1 — p c showing the re- 
manent magnetization with the usual order-parameter 
behavior [5]- [7] and prominent peak in temperature de- 
pendence of magnetic susceptibility which appears in 
low fields as a result of dilution and becomes gradually 
smeared in higher fields [EMU]- 

To get some notion of the thermodynamics of possi- 
ble DAFF ferrimagnetism we consider here the nearest- 
neighbor AF on the simplest hierarchical lattice, imitat- 
ing the percolation cluster with fractal dimension d = 2 
and 77 = 1/3. It may describe qualitatively dilute AF, 
large planar aggregates of AF particles or disordered AF 
thin films. The model also exhibits a number of field- 
induced ground state transitions marked by the mag- 
netization jumps which are characteristic for dilute AF 

[23, El- 



LOW-FIELD THERMODYNAMICS OF 
HIERARCHICAL ANTIFERROMAGNET 

We consider the short-range Ising AF on the sim- 
plest "diamond" hierarchical lattice [T2]. Its building 
process is shown in FigJ5J On the 77-th level of hier- 
archy the lattice has N n sites, N n = |(4™ + 2), see 
Ref. 13 . The coordination numbers of the sites are 
the powers of 2: z = 2,4,8, .... At all levels of the hi- 
erarchy the lattice is bipartite and for n > the sites 
with coordination number z = 2 constitutes the sub- 
lattice A (open circles in FigJ5]) while the others be- 
long to the sublattice B (filled circles), N A , n ^ Nb, u - 
At the Ti-th level N A ,n = 2 • 4"" 1 for n > [15]. so 
Vn = (N n - N A>n )/N A>n = (1 + 2 • 4 1 -") /3. We are 




n=0 n=1 n=2 

Fig. 2. Construction of hierarchical lattice. It is bipar- 
tite at all levels. Different circles designate its partition- 
ing, open circles sublattice A, filled circles sublattice 
B, N A > N B . 

interested in the thermodynamic limit of infinite levels 
of hierarchy. In this limit 77 = 1/3 and fractal dimension 
d = 2 [T5J. 
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For the Ising spins Si = ±1 placed on the sites of 
this lattice we consider the AF Hamiltonian 

n = J Yl SiSj-HAj^Si-HB^Sj (1) 



<ieA,jeB> 



ieA 



j£B 



where < i € A, j € B > means the summation over 
nearest neighbors and different fields for the sublattices 
are introduced. This allows to find the average mag- 
netizations of each sublattice and the order parameter 
for the transition. Homogeneous field corresponds to 
H A = H B = H. 

The usual way to get the partition function of the 
model is through the recursion relations for partial par- 
tition functions at different levels of hierarchy Z n (S, S') 
having fixed values of the outmost left and right spins 
S and S' [El. These relations read 



Z n +i(S, S') 



]T Z n {S,S 1 )e h ^Z n {S 1 ,S f ) 



.Si=±l 



(2) 

h n = H n /T, H = H A and H n = H B> n > 1. The 
initial condition for them is 



K = J/T. (3) 



Z Q (S,S') = e- KSS ', 
Using the representation 

Z n (S, S') = cxp ~ [C n + u n SS' + K - h n ) (S + S')} 

we get from Eqs.(2,3) 

u = -2K, v = h A , C = 
u n +i = 2 In cosh u n + In (l — tanh 2 u„ tanh 2 v„) , 
w„ + i = 2v n + 2tanh _1 (tanh u n tanh v n ) — 2h n + h n +i, 
C n+ i = 4C„ + u n+1 + 4 In (2 cosh v n ) . (4) 

The last of Eqs.(4) gives for n > 

n-l 

C n = U n - 4™U + ^ 4 "~' l U l + ln ( 2 cosnu i)] 
1=0 

so the n-th level free energy per spin is 

F n = -^-lnJ2Zn(S,S')e h °( s+s ') = 

" S,S' 

3 o 

- -TJ2 4"' N + In (2 coshv,)] -~J + (1/N n ) 



1=0 



(5) 



At H n = the model has phase transition at K 
K c « 0.609 being the solution to the equation K c 



lncosh2i4r c . u c = 2K C is the stationary point of the 
zero-field equations, u n +i = 2 In cosh u„, v n = 0. In the 
paramagnetic phase at K < K c u n — > for n — > oo, 
while in the ordered phase at K > K c u n — > oo. Ac- 
cording to above considerations the ground states of the 
model have magnetizations ± (1/2), so we may expect 
that the ordered phase is fcrrimagnetic. To show this 
we consider Eqs.(4) at 

0< (T c - T)/T c = r « 1, |Ii„|«t. (6) 

In this and v n can be found approximately in 

the three regions of n: 

1) 1 < n < A, u n - u c < u c , \v n \ < 1 

U n ~ U c + (Ul ~ U c ) , (7) 

h B + {2 + k)^ 1 nh - 
v„ ~ — , h = h B - (k + 1) hA (8) 

1 + K 

ui = 21ncosh2if, k = 2tanh2iT c « 1.68, 
The value of A is defined by 

ma = 2u c , k~ a = — « r (9) 

while |i> n | <C 1 requires 

\h\(2 + n) X = |^| T -M2+«)/i»ic ^ i ( 10 ) 

2) A < n ^ /i, 1 « u.„ « |w n | 

n-l 

u n ~ 2 n - x+1 u c - Y 2 n " fc ln2coshw fc , (11) 

k=\ 



Vn^^+r- X [v X -^), (12) 



^A 



Kb 

:-! V ~3 

(1 + «) ~ (1 + k) (2 + k) r ln (2+«)/ln« 



{2 + n) X - 1 nh 



(13) 



The value of fi is defined by the equation u^ = 
As u M w 2^- A+1 u c - £ 2^ l - fe ln2cosh4' £ - A v A w 

fc=A 

2^ A 2 (u c - In 2) + 4^ A t;A, v n rj 4^" a i;a, we get 
„ _ (u^n^, _ _ 

3)/i > n,|« n | > l,u„ w 0,v n 2 n "' 1 u /1 . (15) 

Note that in the sums we consider the large numbers A 
and u as integers neglecting its fractional parts. 

Using the above approximations for u n and v n we 
can find from Eq.(5) free energy in the thermodynamic 
limit near the transition point in a small field (cf. 
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Eqs.(6), (10)). Thus, dividing the sum in (5) in three 
parts and n — term, 



we get 



AF 

--- = In 2 cosh /i A + £ A + S Aai + S M , 



A 

S A = ^[4-" ( Uc + ln2)+ Uc r(|) 

n=l 

+ 5i2(l + «) 2 (2 + «) 2 ^ + 2) 



2)1 



3(^2 + ^) + ^;^- — 



+ - 



2u c + In 2 

4 — AC 

K/l 2 / K.\ 2A 



On)' 



2(4 + k)(1 + k) 
E v = ^ (2-"- A+1 Mc + 4-"ln2coshw„)- 

n=A+l 
M n— 1 

^ 4-™^2 n - fe ln2coshi; fc 

n= A +l fe= A 

= ^ (2-"- A+1 u c + 4-"ln2cosht; n )- 

n=A+l 

(4" fc - 2- k -f) In 2 cosh u fe 

k=\ 

- 4~ x 2u c - 4"^ + 4~ M In 2 cosh^ - 4~ A In 2 cosh v x 
:4- x (2u c -ln2) 



E 4-" K| « 4-^K| = 4"> A |. 

n— 

Here we used |u A | < 1, > 1, Eqs.(7, 8, 11, 12, 
15) and the relation following from Eq.(ll), 

^2- fe -^ln2coshw fe = 2- x -" +1 u c -4"^. 

fe=A 

Finally we have from Eqs.(5, 9, 13, 14) 



F/T c ss - X c / 2 - In 2 + rs c - aT 2 ~ Q - br' 3 



-ct^ 2 , 
(16) 

s c = In 2 - 2K c tp± « 0.34, a = 2^^^ - In 2 w 0.17, 

6 = 4(2+kK1 + «) ~ °' 13 ' C = 8(4+«0a + «) a ~ °' 015 - 



„ In 4 „ „„ n In4-In(2 + K) _ 

a = 2 - - — w -0.67, /3 = — - w 0.16, 

In k In k 

7 ^ 21 ^ 2 + ^- lD4 .2.35 (17) 
ln« 



In homogeneous field h = —nh (cf. Eq.(8)) so F in 
Eq.(16) has the standard scaling form of a ferromag- 
net with spontaneous magnetization m ~ t 13 and diver- 
gent susceptibility \ ~ t~ 7 . This expression is valid at 
< r < 1, |/i| < t^+t, cf. Eq.(10). Scaling indices 
(17) obey the usual relation a + 2/3 + 7 = 2. Negative a 
means that specific heat is finite at the transition point 
and has a cusp at T c . Note also that s c is the entropy at 
the transition point. So this AF system looks like gen- 
uine ferromagnet, even featuring the absence (smearing) 
of transition in a finite field. The last is evident as the 
nontrivial stationary point of finite-field recursion rela- 
tions (4) cannot be reached from any initial conditions. 

Yet the dependence of F on h from Eq.(8) shows that 
true order parameter for the transition is a linear com- 
bination of Ma — E &i an d Mb — E &i conjugate 

with H = Hb — (k + 1) Ha- To distinguish the order 
parameter in the Hamiltonian (1) we perform a coordi- 
nate rotation in 2d space of vectors H = {Ha, Hb) and 
M = (M A , M B ) to bring the term -MH in (1) to the 
form -MH = -MH - M'H' where 



^ = M B - (k+1)M a m , _ (k + 1) M B + M A 



i + (K + iy 

H' = {k+1)H b +H a 



Thus M is the order parameter while M' and H' are 
non-critical variables. Hence, (M') = at H — > so the 
spontaneous magnetic moments of the sublattices obey 
the relation (Ma) = — (k + 1) (Mb)- Then for the spon- 
taneous magnetizations m„ — (M v ) /N v (u = A, B) 
we have 

m,A = —i](k + 1) rag w — 0.9to_b 

This differs from the ground state relation vtia = — tub- 
We may suggest that this is a consequence of critical 
fluctuations diminishing wia more strongly than tub 
as all sites of sublattice A have the lowest coordina- 
tion number za = 2. To some extent this effect would 
be present in all dilute AF on imbalanced bipartite 
graphs since the sublattice A with larger amount of 
spins would necessary have lower average coordination 
number za — C/Na < zb — C/Nb- Here C is the num- 
ber of bonds and we used the fact that all bonds are of 
A-B type. Thus za = i]Zb so the lower 77 the stronger 
can be the fluctuation-induced disbalance between vtia 
and in b near T c . At r/ = 1 this effect vanishes so its 
observation in neutron-diffraction experiments can cer- 
tify the onset of imbalance in the magnetic percolation 
cluster. 
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GROUND STATE TRANSITIONS 

Here we assume Ha = Hb = H. At T = wc define 
u n = lim Tu n , v n — lim Tv n , E = lim F to obtain 

T^O T^O T->0 

from Eqs.(4, 5) 

Uq = — 2J €'o = H 
u n+1 = 2 (|«„| - \v n \)ti(\u n \ - \v„\) , (18) 
v n+ i = 2v n + 2min (|u„| , |tT„|) sgn (u n v n ) - H, (19) 

OO 

-(4/3)E = H + Y / 4- n (u n + \in\) 

n=l 

i9 in Eq.(18) is Heaviside's step function. Solving 
Eqs.(18, 19), we get 



dE _ 3 
dH ~ 4 



(4/3)E = H + J2^ n \H - 2 n+1 J 

n=l 

oo 

I + ^V ,l sgn(tf-2" +1 j) 



So at i? ^ 2 fc J 

m = 2 _1 t? (2 J — ff) + (l — 2 • 4~ r ) §(H — 2 J) , 

where r = [log 2 (H/J)] is an integer part of log 2 (H/J). 
At i? r = 2 r J, r > 2, we have m r = 1 - 6 ■ 4~ r . Field 

m 1 



0.6 



0.4 



*0 4 8 12 16 H/J 
Fig. 3. Field dependence of the ground state magnetization. 

dependence of the ground state magnetization is shown 
in Fig{3j Due to the imbalance ratio rj = 1/3 the system 
has spontaneous magnetization m = 1/2 at H — > +0. 
The jumps at H r = 2 r J result from the flipping along 
the field of single spins in sublattice B having the co- 
ordination number 2 r . Note that in dilute crystalline 
lattices there are many more jumps appearing at ratio- 
nal values of H/J where flipping of the whole clusters 
takes place [TU]. Such jumps were observed in low-T 
experiments in Fe p Zn\- p F2 [TT]. The ground states 
of DAFF on cubic and bec lattices are studied numer- 
ically for all p in Ref. [T3]. Unfortunately, the data 
on the sample's and percolation cluster magnetization 
are totally absent in this paper so we cannot decide if 
there is the imbalance in real 3d percolation clusters at 
p c < p < 1 -p c . 



DISCUSSION AND CONCLUSIONS 



The decades of experimental investigations of dilute 
Ising AF have shown that DAFF - RFIM correspon- 
dence works reasonably well at low dilution [TS] . Mean- 
while the field-induced rounding of the transition ap- 
pears at lower p which is impossible in the case of the 
AF ordered phase. One explanation assumes that this is 
nonequilibrium effect due to the pinning of AF domain 
walls by the vacancies which results in very slow relax- 
ation to the equilibrium AF structure [TB] . The numer- 
ical simulations demonstrate that this effect is strongly 
enhanced with the appearance of vacancies' percolation 

Yet the experiments show that not just a slowing 
down but true nonergodic irreversibility with multiple 
metastable states sets in at intermediate p [18], so one 
may suggest that AF transition transforms into a spin- 
glass one [TSJ,[T5]- This is quite probable scenario at 
p close to p c where weak next-nearest interactions may 
cause a frustration in a sparse percolation cluster. How- 
ever, the field-induced rounding and the nonergodic ir- 
reversibility sets in right below p = 1 — p c in real dilute 
AF gD]. The observation of Ref. [10] that AF transition 
breaks up with the appearance of vacancies' percolation, 
along with the ferromagnetic anomalies found in Refs. 
[5] [9] at p c < p < 1 — pc, strongly supports the on- 
set of imbalance in magnetic percolation cluster mak- 
ing the transition ferrimagnetic in this concentration 
range. This explains the field-induced rounding, while 
the weak dipole-dipole interactions would give rise to 
a number of metastable magnetized domain states thus 
explaining the appearing nonergodicity. The expected 
ferrimagnetic phase at intermediate dilution would dif- 
fer macroscopically from the spin-glass one by a small 
scaling index of remanent magnetization j3 < 0.5 and 
by the divergent magnetic susceptibility. The present 
model also predicts rather high index 7, large negative 
a and disbalance in the sublattice magnetizations near 
T c . But to observe these features of ferrimagnetic tran- 
sition the measurements in ultra-low fields (same as in 
Refs. [S]-[5]) are needed to avoid the smearing of the 
transition. 

Here we can only give qualitative arguments in favor 
of ferrimagnetism of real dilute AF, explaining the re- 
sults of Refs. [5] -[S], [H], [20] ■ Now we do not know how 
to rigorously prove that percolation cluster on dilute 3d 
lattices has r\ v < 1 at p c < p < 1 — p c . Yet this hy- 
pothesis can be verified by numerical simulations. This 
and further low-field experiments may help to elucidate 
the nature of transition in nearest-neighbor dilute AF. 
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